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ABSTRACT 


Time profiles of gamma ray bursts (GRBs) are ex¬ 
tremely diverse in their durations, morphologies, and 
complexity. Nevertheless, the average peak-aligned 
profile of all bursts detected by BATSE with suffi¬ 
cient data quality has a simple “stretched” exponen¬ 
tial shape, F oc exp[—where t is the time 
measured from the time for the peak flux, Up, of the 
event, and to is a time constant. We study the be¬ 
haviour of to of both the post-peak and the pre-peak 
slopes of the average time profile as a function of the 
peak brightness range of the burst sample. We found 
that the post-peak slope shows time dilation when 
comparing bright and dim bursts, while the pre-peak 
slope hardly changes. Thus dimmer bursts have a 
different shape - they are more asymmetric. This 
shape-brightness correlation is observed at a 99.6% 
confidence level. Such a correlation has a natural ex¬ 
planation within the pulse avalanche model, which 
is briefly described. Complex events, consisting of 
many pulses are more symmetric and are intrinsically 
brighter. Bursts consisting of one or a few pulses are 
intrinsically weaker and more asymmetric. For such 
a correlation to be observable requires that the lumi¬ 
nosity distance distribution of GRBs to be different 
from a power-law. 

Keywords: Gamma-ray bursts, Methods: Data anal¬ 
ysis. 


1. INTRODUCTION 


Stern (1996) found that the average peak-aligned 
time profile of GRBs (the procedure of peak- 
alignment was pioneered by Mitrofanov et al. 
1994, 1995) in the BATSE-2 catalog has a simple 
“stretched” exponential shape, 


< F/Fp >= exp[—(t/to)^^^], where t is the time since 
the peak flux, Fp, of the event, and to is a constant 
ranging from 0.3 s for strong bursts to ^ 1 s for 
dim bursts. This dependence of to on brightness can 
be interpreted as a cosmolo Efical time dilation (e.g., 
Paczyhski 1992 ; Piran 19^ ). 


Such a simple average time profile is remarkable con¬ 


sidering the diverse and chaotic behavior of the in¬ 
dividual time profiles of GRBs. On the other hand, 
the simple shape of the average time profile gives an 
excellent opportunity to study effects such as time 
dilation. 

Here, we study the two slopes of the average time 
profile for a larger sample of GRBs and with a more 
accurate treatment of the background than was done 
in Stern (1996). Another advantage we now have is 
the access to the pulse avalanche model developed 
by Stern & Svensson (1996) which successfully de¬ 
scribes many statistical properties of GRBs including 
the stretched exponential shape of the profile and, of 
particular importance for the present work, the rms 
variance of individual time profiles. This means that 
we can rely on this model when estimating the er¬ 
rors of stretched exponential fits, which in turn gives 
us reliable estimates of the significance levels of the 
observed effects. 


2. DATA PROGESSING 


This work is based on data obtained from the pub¬ 
licly available BATSE database at Goddard Space 
Flight Genter. Our sample includes bursts up to trig¬ 
ger number 3745. We used the 0.064 s and 1.024 s 
time resolution data from the Large Area Detectors 
(LAD). All time profiles were constructed with 64 ms 
time resolution together with pre and post-burst ex¬ 
tensions of 1024 ms time resolution. All background 
fits were done with the 1024 ms data as they cover a 
wider time interval including the pre-trigger history. 
The time profiles were studied using count rates in 
all 4 lad’s energy channels. 

The procedure of background fitting included: 

- A visual examination of all bursts including both 
64ms and 1 s resolution data. 

- All doubtful peaks and count rate variations with 

exceeding that of Poisson noise were analyzed in 
order to see whether they came from the same di¬ 
rection as the main peak of the burst. This was ac¬ 
complished by comparing the direction of the eight- 
component vector consisting of the from the eight 









LADs for the feature with the direction of the corre¬ 
sponding vector for the main peak. 


- Discarding all events with highly variable back¬ 
ground as well as those where we were unable to con¬ 
fidently extract sufficiently long intervals of pre-peak 
and post-peak histories. 

- Possible use of widely separated fitting windows to 
avoid losses of weak GRB signals. 

When sorting bursts into brightness groups we used 
peak fluxes for 64 ms time resolution from the 
BATSE-3 catalog. 

The errors were calculated using the pulse avalanche 
model. What we need are the statistical errors of the 
stretched exponential fits to the data. To extract this 
error directly from the data is very difficult because 
of strong correlations along the time profile. Instead, 
we simulated many samples of N events, determined 
the average time profile for each sample, and then 
made a stretched exponential fit to each simulated 
average profile: F{t) = /3exp[—where to 
and /3 are fitting parameters. The same procedure 
was used as for real data. Finally, the rms errors of 
the fitting parameters was calculated. 

We have two slopes of the profile - the pre-peak (ris¬ 
ing) slope and post-peak (decaying) slope. We fitted 
them simultaneously with different to, denoted here 
as tr and td, and with a common (3. Fitting the sim¬ 
ulated samples gave the following relative standard 
deviations: 

(T{tr4)ltr,d = 0.20v/TW^ 


a(tr -f td)/{tr + td) — O.IQa^/IOO/A^ 


<j{tdltr)litdltr) = o.iayiooTiv 

Note, that the relative accuracy for the sum of the 
two time constants is close to that for one time con¬ 
stant, while the accuracy of the asymmetry ratio, 
td/tr, is markedly better. This is a consequence 
of a strong correlation between the two slopes - a 
circumstance that favors the measurement of shape 
vs. brightness correlations and that complicates the 
measurement of the time dilation effect. The errors 
are robust against variations of the parameters of the 
pulse avalanche model as long as the model gives ap¬ 
proximately the correct stretched exponential aver¬ 
age profile. The procedure for the data analysis will 
be described in greater detail in Stern et al. (1997). 


3. THF AVFRAGF TIME PROFILE FOR 
DIFFERENT BRIGHTNESS GROUPS 


The average peak-aligned profiles for three brightness 
groups are shown in Figure 1. Both the rising and 
and the decaying profiles are well shaped stretched 
exponentials for the bright and the medium group, 
while both profiles are quite deformed for the weak¬ 
est group. The rising (pre-peak) slope is steeper for 
all brightness groups, but the asymmetry is increas¬ 
ing when going from the brightest to the dimmest 
group. The results of our stretched exponential fits 
are summarized in Table 1 and in Figure 2. 



Figure 1: Average peak-aligned time profiles for three 
brightness groups: 1.) Fp > 5 ph cm~‘^ 157 GRBs 

(upper eurves); 2.) 0.7 < Fp < 5, 630 GRBs (middle 
eurves); and 3.) Fp < 0.7, 116 GRBs (lower curves). 
Middle and lower curves are shifted downwards for clar¬ 
ity. Solid and dotted curves represent the average post¬ 
peak and pre-peak time profiles, respectively. Straight 
lines show the best linear fit to the post-peak history of 
the medium brightness group. 



Figure 2: Time constants, tr and td, vs. peak photon flux, 
Fp, in 64 ms time resolution. Lower and upper crosses 
represent tr and td for the pre-peak (rising) and post-peak 
(decaying) average time profile, respectively. Error bars 
of the time constants correspond to la. Error bars in 
photon flux represent the width of the brightness groups. 











# 

Peak flux 

N 

td 

tr 

tr + td 

td/tr 

1 

12.5 - 200 

64 

0.42± 0.10 

0.36±0.09 

0.78 ± 0.18 

1.17 ±0.19 

2 

3 - 12.5 

193 

0.60± 0.08 

0.40±0.06 

1.00 ± 0.13 

1.50 ±0.11 

3 

1.75 - 3 

159 

0.75± 0.11 

0.37±0.06 

1.12 ± 0.16 

2.02 ±0.21 

4 

1 - 1.75 

241 

0.72± 0.09 

0.46±0.06 

1.18 ± 0.14 

1.56 ±0.13 

5 

.7 - 1 

139 

0.85± 0.14 

0.48±0.08 

1.33 ± 0.21 

1.89 ±0.18 

6 

0 - .7 

116 

0.80± 0.14 

0.37±0.07 

1.17 ± 0.20 

2.16 ±0.26 

7 

7.5 - 200 

111 

0.47± 0.08 

0.36±0.07 

0.83± 0.14 

1.30 ±0.16 

8 

5 - 200 

157 

0.50± 0.07 

0.35±0.06 

0.85 ± 0.12 

1.56 ±0.13 

9 

.75 - 2.5 

463 

0.78± 0.07 

0.45±0.06 

1.22 ± 0.11 

1.73 ±0.10 


Table 1: Time constants [s] of the stretched exponential fit to the averaged pre-peak (U) and post-peak (td) 
profiles. Peak flux [ph cm“^ s“^] is taken from BATSE database. N is the number of bursts in the given 
brightness interval. 


Fits of all time profiles have a good except that for 
the weakest sample (# 6 ). Its profile is apparently 
deformed. This deformation is partially associated 
with the low trigger efficiency of this sample and with 
effects of Poisson noise. Nevertheless these biases are 
insufficient to explain this deformation and a real de¬ 
formation effect could be present in the shape of the 
weakest sample. An analysis of this effect is beyond 
the scope of this poster paper. 

The strongest effect in the remaining samples is the 
time dilation of the decaying (post-peak) slope of the 
time profile. Comparing td for samples 1 and 9 we 
find the time dilation to be a factor 1 . 86^0 5 ?) where 
the 90% confidence interval is given. Comparing sam¬ 
ples 7 and 9 gives a factor I.OGIq;®®. Varying the 
lower flux limit of the brightest sample from Fp = 
5 ph cm“^ s“^ and upwards does not significantly 
change the magnitude of the time dilation effect. A 
careful estimate of the significance level of the time 
dilation using model simulations gives 0.985. 

The rising (pre-peak) slope is, however, surprisingly 
stable. The variations of tr do not exceed statistical 
errors. This leads to an increasing asymmetry ratio, 
td/tr, for the weaker samples. The effect is significant 
comparing samples 1 and 9: (td/tr)dim/{td/tr)hright 
= l-48to;32 (90% confidence interval). The proba¬ 
bility of sampling such ratios if the asymmetry was 
constant is 4 • 10“^. 

The time dilation of F + td has a smaller, but still ac¬ 
ceptable significance when comparing samples 7 and 
9. The “dim/bright” ratio becomes 1 . 47 tQ; 3 |. The 
probability for zero time dilation is 0 . 01 . 


4. INTERPRETATION OF CORRELATIONS 
USING THE PULSE AVALANCHE MODEL 


The pulse avalanche model is based on the assump¬ 
tion that the well-shaped stretched exponential time 
profile results from a simple stochastic process re¬ 
sponsible for the generation of time histories of 
GRBs. Then the diversity of CRB’s time histories 
arise as different random realizations of the same 
stochastic process at approximately the same param¬ 
eters. The important requirement is a near-criticality 


of the process - then it provides a large variety of in¬ 
dividual bursts behaviours. 

This idea was implemented as a near-critical chain 
reaction of events (of still unknown nature, it could, 
e.g., be reconnections of turbulent magnetic field), 
where each event is associated with one pulse of 
gamma-ray emission. Then, in a near critical regime, 
one spontaneous pulse can give rise to a long cas¬ 
cade of secondary pulses piling up into a complex 
chaotic event, or, depending on chance, no further 
pulses may result and we will instead see a simple 
single pulse event. 

With a proper scaling of time delays between 
pulses and a proper spectrum of spontaneous pulses 
(flicker noise), the model successfully reproduces the 
stretched exponential profile (see Figure 3) as well 
as the autocorrelation function. There is also quali¬ 
tative agreement with the duration distribution, the 
ratio of the number of simple and complex bursts is 
reproduced, and the model even produces the visual 
impression of real bursts. For more details see Stern 
& Svensson (1996). 

Just two simple assumptions of those implemented 
in the pulse avalanche model are needed to demon¬ 
strate that all correlations described above are very 
natural. GRBs consist of a number of pulses of differ¬ 
ent durations but similar shapes. Let all these pulses 
have locally independent sources of energy. Then if 
two pulses coincide in time, their amplitudes sum up. 
In a complex event, hundreds of pulses are piling up, 
increasing the peak brightness by up to an order of 
magnitude. Then, if a pulse is a kind of standard can¬ 
dle, simple events are intrinsically weak and complex 
events are intrinsically bright. 

At the same time, simple events are asymmetric just 
because a single pulse is asymmetric (e.g., Norris et 
al. 1996) with a sharp rise and a slower decay. This 
asymmetry is washed out in complex events where 
the position of the highest peak is more or less ran¬ 
dom among many overlapping pulses. 

To demonstrate this we simulated a large sample of 
“bursts” using model parameters which gave an ap¬ 
proximate agreement with the stretched exponential 
average time profile for the sample of all real bursts. 
The results are summarized in Table 2. We see that 
the asymmetry ratio, td/tr, increases with decreasing 
peak flux. 



Peak flux 

td 

tr 

tr + td 

td/tr 

U - oo 

EI^ 

0.41 

"TUT 

1.61 

3-00 

0.75 

0.68 

1.43 

1.10 

0.8 - 3 

0.72 

0.50 

1.22 

1.44 

0 - 0.8 

0.50 

0.19 

0.69 

2.63 


Table 2: Time constants of the simulated bursts in 
different “intrinsic brightness” intervals. The ampli¬ 
tude of the single pulse is sampled uniformly from 
the interval [0,1]. 
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Figure 3: Comparison between observed and simulated 
average post-peak time profiles of GRBs. Thin-line his¬ 
togram labelled < F/Fp > shows the average peak-aligned 
post-peak time profile for the 598 useful BATSE-3 events 
as the fractional flux, < F/Fp > vs. where Fp 

is the peak flux and t = 0 - 150 s is the time since 
the strongest peak. Thiek curve labelled < F/Fp > 
is the average peak-aligned time profile for 5000 simu¬ 
lated time profiles. The set of curves labelled Tl-ms O'l'e 
the rms deviations of individual peak-aligned time pro¬ 
files, Frms = [< {F/FpY > - < F/Fp for 

both real (thin line histogram, 598 events) and simulated 
(thick curve, 5000 events) time profiles. Further details 
are given by Stern (1996) and Stern & Svensson (1996). 


5. CONCLUSIONS 


Besides the time dilation effect, observed in many 
previous works (e.g., Norris et al. 1994) we also see 
a dependence of the profile asymmetry on brightness 
and this effect is of the same order of magnitude as 
the time dilation itself. 

A correlation of such a kind cannot be due to spectral 
redshifts. We hnd that strong events have slightly 
smaller asymmetry in the higher energy bands, (LAD 
channels 4 and 3) than in the lower energy bands 
(LAD channels 1 and 2). Redshifting the softer part 
of the spectrum below the detector threshold would 
then give rise to more symmetric, rather than asym¬ 
metric profiles. Details will be published elsewhere. 

The effect of trigger efficiency is negligible for the 
brightness groups considered here. Also one can 


hardly find an evolutionary factor that would change 
the asymmetry. One can suggest that there exist 
two separate classes of GRBs with different degree of 
asymmetry, which are differently distributed in space 
(separate classes of long bursts are required as short 
bursts do not contribute to the asymmetry). But 
such a suggestion seems too arbitrary, too radical, 
and unnecessary as there exists a much simpler ex¬ 
planation. 

The simplest explanation is that the observed corre¬ 
lation is a consequence of an intrinsic correlation be¬ 
tween shape and brightness as described above. The 
necessary condition for such a correlation to be ob¬ 
servable is a significant deviation from a power law for 
the CRB distribution over luminosity distance. This 
would allow intrinsically strong events to dominate 
in the brightest observational range. The observed 
log N - log P distribution is actually curved (Mee- 
gan et al. 1996) and this is natural if the distance 
distribution covers both Euclidean and z ~ 1 regions 
which have different luminosity distance scalings. 

Maybe our detection of a strong shape - brightness 
correlation imposes a stronger constraint on the cur¬ 
vature of the true radial distribution of GRBs than 
what follows from the observed log N - log P dis¬ 
tribution. Detailed studies are, however, required to 
formulate this intuitive conclusion at a quantitative 
level. 

As one kind of correlation has been observed, other 
kinds of intrinsic correlations may also be observable 
and this causes a problem for the cosmological in¬ 
terpretation of the time dilation effect. This prob¬ 
lem, considered by Brainerd (1994), arises from un¬ 
avoidable correlations between peak luminosity and 
time scales caused by different bulk Lorentz factor in 
the sources of different bursts. This effect can mimic 
both cosmological time dilation and spectral redshift. 

However, if our interpretation using the pulse 
avalanche model is valid we must conclude that the 
real time dilation is larger than that obtained from 
Table 1. Actually, within the pulse avalanche frame¬ 
work, intrinsically weak events are not only more 
asymmetric, but they are also narrower (see column 
tr + td in Table 2). This is a new correction that 
increases the real time dilation more than the cor¬ 
rection arising from spectral redshift (see Norris et 
al. 1994). The corrected time dilation could exceed 
a factor 2 and it could be caused by different effects, 
including the cosmological one. Unfortunately, the 
task of extracting the cosmological component from 
the total time dilation seems extremely difficult. 
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